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Abstract 

This study is devoted to the two major hydro-ecosystems of 
the Carpathian region in Ukraine, the Dniester and Danube 
Rivers. For the first time, the patterns of changes in the 
ecosystems’ quality parameters were established by means of 
developing functional dependences among the ecological 
standard values of the composite quality indicators of natural 
waters and the river length and terrain altitude. Using the 
statistical processing of the data, the quality monitoring of the 
upper reaches of Dniester and Danube ecosystems in the 
Ukrainian Carpathians was conducted for the period from 2001 
to 2019. The data on the Composite Potential Quality 
Indicators standard values was gathered for the heights of the 
Carpathian region: at 50 m interval for the altitude up to 1,000 
m, and at 100 m interval for more than 1,000 m altitude. The 
acquired dependencies can be used as the territorial 
background standard values of the ecological indicators of 
basin ecosystems. It will help to scientifically substantiate the 
ecologically safe values of the anthropogenic pressure. 
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Introduction 


According to the Ukrainian environmental legislation! (Law of 2017), environmental quality assessment, 
including water resources, is carried out to set the maximum permissible levels of various environmental 
indicators. It guarantees the environmental safety to the population and ensures natural resource 
management, genetic conservation and sustainable development. 


The Ukrainian Carpathians comprise the basins of four major rivers (Mandryk et al., 2017) that are grouped 
into two regional hydro-ecosystems (Figure |): Dniester (right-bank tributaries 1n the territories of Lviv and 
Ivano-Frankivsk oblasts) and the Danube rivers. Later includes the Prut river within the territories of Ivano- 
Frankivsk and Chernivtsi oblasts, together with the Siret tributary in Chernivtsi oblast and the Tysa tributary 
in Transcarpathian oblast. These waterways are of international importance, as their basins are located 
within several countries. The Prut and the Siret rivers flow through Ukraine, Romania and Moldova; the 
Tysa river crosses through Romania, Hungary and Slovakia; the Dnister river flows through Moldova 
(Kinash et al., 2019). For the reason of cross-boundary water flows, this investigation and modeling of the 
upstream quality of international water flows becomes an important task (Odnorih et al., 2020). 


The watershed of the Ukrainian Carpathians covers about 28,000 watercourses, most of which are small 
ones with a total length of more than 50,000 km. The river network density averages 1-1.5 km per km? 
(Prykhodko et al., 2020). Due to the mountainous terrain of the Carpathian ranges, the basin areas are small. 
Small rivers (up to 10 km long) dominate and account for almost 98% of all watercourses (Khilchevsky1, 
Kurylo and Sherstyuk, 2018). There are 457 rivers over 10 km in length. In the Prut river basin, all rivers, 
except for the main river, are classified as small (Karpinski et al. 2018). There are seven rivers having 100 
to 300 km length in the basins of the Dniester, Tysa and Siret rivers. The total area of river surface basins 
in the Carpathians exceeds 37,500 km? (Mandryk et al., 2020). 


The purpose of the study is to establish the background values of the qualitative parameters of the hydro- 
ecosystems of the Ukrainian Carpathians and to determine the ecological standards of the components of 
natural waters within the study area. To achieve this, the general and individual patterns of their spatial 
distribution based on statistical processing of data from the results of hydro-chemical observations has been 
studied, using data from the state network for monitoring the quality of water bodies for the period 2001-2019. 


Methodology 


To undertake the study, following research methods were used. First was an analysis of environmental 
information on hydro-ecosystems and its processing; and second was the statistical and mathematical 
analysis using Microsoft Excel and TableCurve 2D software (Stevens, Springer and Ledbetter, 2011). 


One of the complex criteria used to assess water quality in Ukraine is the Water Pollution Index (WPI). The 
WPI is a composite quantitative assessment of water quality by a set of key indicators and types of water 
use (Singh, Yadav and Yadava, 2016). This integral index is calculated by the following formula: 


WPI = 3(C/ MPC)/n, (1) 


where MPC is maximum permissible concentration (value) of the indicator; 
C is actual concentration (value) of the indicator; and 
nis number of indicators. 


1 


https://zakon.rada.gov.ua/laws/main/en/tl 
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Legend: 


[Sex] . hydro-ecosystem of the Dnister River 
[Se] - hydro-ecosystem of the Prut River 

[Nex | - hydro-ecosystem of the Tysa River 
[*Pe| - monitoring network of the Dnister River 
_@®5| - monitoring network of the Prut River 
[+T @| - monitoring network of the Tysa River 
. the Carpathian region 
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Figure 1: State monitoring network of water bodies of the Carpathian region of Ukraine 


For surface waters, the number of indicators used to calculate the WPI should not be less than 5 regardless 
of whether they exceed the MPC or not. But they must include the dissolved oxygen (DO) and biological 
oxygen demand for 5 days (BODs) (Shmandly et al., 2017). The advantages of this index are the ease of 
calculation, data accessibility and availability. That is why it is most often used in the State Water Resources 
Agency of Ukraine (the central executive body, which implements the state policy in the field of water 
management development) to assess river pollution. 


Despite popular usage of WPI, the WPI does not sufficiently assess clean rivers. In the classification, there 
is quality class 1 (very clean water) having WPI < 0.2, and class 2 (clean water) having WPI from 0.2 to 1 
(Government of Ukraine, 2018). It is believed that the Water Pollution Index is not flexible and sufficient 
to assess the degree of self-purification of rivers, the dynamics of changes in the ecological potential 
depending on the height and width of the terrain, water content of rivers, dynamics of changes in the 
anthropogenic load, and so on. 


Considering the above-explained limitations of WPI, the Composite Potential Quality Indicator (CPQI) 
(Arkhypova et al., 2019) was used to determine natural regularities in this article. The CPQI is obtained by 
adding the coefficients of indicator reserves (relative reserve capacity) if the actual value is less than the 
permissible value, and by subtracting if the actual value is greater than the allowablepermissible value. The 
product is divided by the number of values used. The excess of the permissible value on the actual value 
shows the coefficient of stock (potential, reserve capacity) of the hydroecosystem. Moreover, the excess of 
the actual value over the permissible value shows the stock deficit ratio (relative reserve shortage). 


The outcome is divided by the number of indicators used. 
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where QS; is the water quality standard for the 7” indicator—limit values of water status indicators that 
meet the requirements of different types of water use; 

C; is actual value of water quality for the i” indicator; 

n is the number of indicators. 


The water samples for physicochemical studies were taken at least 4 times a year (in different hydrological 
seasons) in accordance with the procedure prescribed by the normative document KND 211.1.1.106-2003 
(Zasidko et al., 2019). The following indicators were used to calculate the WPI: NH4* (MPC — 0.5 mg/dm°), 
NOz (MPC — 0.08 mg/dm*), NO (MPC — 40.0 mg/dm?), BODs (MPC — not more than 3 mg/dm7°), dissolved 
oxygen (MPC — not less than 6 mg O/dm’). 


This research used a statistical model based on the mathematical processing of the statistical dataset. The 
data included hydro-chemical analyses of surface waters conducted by the certified laboratories of the 
Dniester Basin Authority and the Analytical Control Department of the Carpathian National Nature Park. 
Each numerical result, corresponding to a specific monitoring section, is an integrated value of over 50 
values, as the samples were taken every year, four times a year. The existing hydro-chemical monitoring 
database was processed by the method of determining the geometric mean: 

F(x) = In x, (3) 


The geometric mean reflects the arithmetic mean of the logarithms of individual values. However, at the 
same time, the effect of the extreme values on the mean is greatly reduced, which is often observed in the 
variation series of monitoring the quality of natural waters (Singh, Yadav and Yadava, 2016). 


TableCurve 2D was used to obtain linear functional dependences and regression equations for statistical 
series of observations. Functional dependence equation, equation coefficients, and calculated statistics 
parameters are represented in the description above the dependences curve. TableCurve 2D uses four 
goodness-of-fit statistics. 


In the following descriptions, SSM is the sum of squares of the mean; SSE — the sum of squares of errors 
(residuals); n is the total number of data values; m is the number of coefficients in the model; DOF is the 
degree of freedom (KreSi¢ and Stevanovic, 2010). 


DOF = n-m, (4) 
where 7° coefficient of determination (r squared): 
r= 1 — SSE/SSM. (5) 
Degree-of-freedom adjusted coefficient of determination: 
DOF r = (1 — SSE*(n-1))ASSM*(DOF-1)). (6) 
Calculated standard errors: 
StdErr = sqrt(SSE/DOF). (7) 
F-statistics: 
F-stat = ((SSM-SSE)Xm-1))AHSSE/DOF). (8) 


The dependence between data sets will be closer if the 7 values approach 1.0 (0 is the total lack of 
correlation), the standard error decreases to zero, and the F-statistics goes to infinity. For the upper reaches 
of the tributaries of the Danube ecosystem, the ecological monitoring database was processed using formula 
(1). The WPI calculations gave long-term values sets of the composite quality indicator, which reveals 
certain patterns of spatial distribution. In hydrometeorological calculations, the concept of standard value 
(as the long-term annual average of the indicator value) is used. For the long-term annual average of WPI 
values (8 monitoring sections at different altitudes), dependence curves and functional linear regression 
equations were obtained that can be the basis for determining the WPI standard value. 
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Results 


The general and individual regularities of spatial distribution of water quality indicators on the basis of 
Statistical processing of data of readings of the hydro-chemical observations were analyzed, using data of 
the state network for monitoring of quality of water objects for the period 2001-2019 (Figure 2). An example 
of calculating the Water Pollution Index is shown in Figure 2, and an example of calculation of the CPQI is 
depicted in Figure 3. 


Research results of the Prut River Water —— Index (2001-2019) 
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Fieure Z, Example of the Water Polluion index calculation with monitoring data processing 
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The pollution level and state of the hydro-ecosystem are assessed based on the value of the Water Pollution 
Index or the Composite Potential Quality Indicator (Table 1). 
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Table 1: Numerical criteria for the aquatic ecosystem state assessment based on the values of composite 
quality indicators 


Water pollution level | WPI value | quality | Aquatic ecosystem state CPOQI value 


Very clean buffer state (ecological balance area) 
Clean 0.2-0.1 
Low pollution 1.0-2.0 
High pollution 4.0-6.0 -3<CPQI<-| 
Very high pollution _| 6.0-10.0 -3 <CPQI<-5 
<-5 


Extreme pollution >10.0 catastrophic state (ecological 
disaster area) 


On the basis of the results shown in Table 1, the analysis was performed for different hydrological seasons 
and monitoring sections. The average annual values were calculated, and the data obtained for the period of 
2001-2019 were summarized. It is observed that the water quality at altitudes above 500 m in the Ukrainian 
Carpathians generally meets the reference conditions. 





. 
2- 
0- 
Medium pollution 
0- 
0- 26: __ 


The long-term monitoring results of surface waters in mountain areas, the beds of which are located at altitudes 


above 500 m, show that the water belongs to the quality classes “clean”, “fairly clean’, “good” or “very good”. 
The discharge of sewage water from mountain settlements diffuses pollution in 99% of the cases. 


For the upper reaches of the tributaries of the Danube ecosystem, the functional dependences among the 
changes in the composite quality indicator of surface water ecosystems and the river length and terrain 
altitude were obtained (Figures 4, 5). The obtained curve (Figure 4) of the dependence between the Water 
Pollution Index standard value and the terrain altitude for the Prut river passes through the centre of the 
cluster of observation points, the curves describing the 95% confidence interval are nearby. Thus, there is a 
close functional dependence between these properties. The coefficient of determination (D=r7=0.82) shows 
that 90% of the variance in Y is due to the variance in X. The actual quality state of the Danube ecosystem 
in Ukraine to the altitude of 700 m is the optimal state, and below this level is the state of adaptation tension. 
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Figure 4: Functional dependence between the change in the Water Pollution Index of the Danube 
ecosystem and the terrain altitude in the Ukrainian Carpathians 
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The significance of the coefficient (77=0.82) was estimated using the F-distribution quantiles tables and the tables 
of coefficient significance. The latter have two entries — the number of degrees of freedom, which is equal to the 
number of observations minus 2; the significance level is accepted at 5%. The critical coefficient value was found 
after determining the significance level and the number of degrees of freedom. In this case, the value of the 
coefficient r°=0.82 was greater than the table critical value r°=0.707 (for sample volume 11 and significance level 
p=0.05) (Mandryk et al., 2020). That is, the null hypothesis of no correlation between the properties was rejected. 
The hypothesis that there is a significant correlation between the properties was accepted. 


Therefore, for the first time, the functional dependence between the change in the long-term annual average 
of the Water Pollution Index in the Danube ecosystem in the Ukrainian Carpathians and the altitude is 
represented. It is expressed by the following significant regression equation: 


WPI = 30.128-H°™, (9) 
where #7 1s the absolute height of the terrain, m 


Despite the continuous diffusion of pollution, which was observed in the Danube ecosystem (even with the 
values of WPI<1), a close relationship between the composite Water Pollution Index and the length of the 
river from its source is evident (Figure 5). 





y = 0.2347e9.0116 
R? = 0.8137 











80 100 
Length from the source, km 





Figure 5: Functional dependence between the quality standard value of the water bodies of the Prut 
ecosystem and their length from the source 


F-test or F-statistics is calculated using the software. The information on the distribution of this value allows 
to check the statistical significance of the regression model based on the value of the coefficient of 
determination. In fact, these tests test the hypothesis that the true coefficient of determination is equal to 
zero. The critical values of F-test were looked up in standard tables. The proposed hypothesis (on the 
existence of close dependence) 1s rejected if the table value of F'stat is greater than the calculated one and it 
is accepted as true if F'stat is less than the calculated one. In all the cases considered, the proposed hypothesis 
was accepted as true if the significance level was assumed to be 1%. In this case, Fstat=56.88 
(calculated)>Fstat=5.99 (table value) (Toms and Lesperance, 2003). Thus, the correlation between the 
properties is not accidental, it is significant. 
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Therefore, as a result, for the first time the linear functional dependence between the change in the Water 
Pollution Index of the Danube ecosystem and the length of the river in the Carpathian National Nature Park 
was found. It is described by the following significant regression equation: 


WPI=0.235-€°09%, (10) 
where L is the length of the watercourse from its source to the monitoring section, km. 


In the conservation areas located in the upper reaches of the Danube ecosystem in the Ukrainian Carpathians, the 
Carpathian National Nature Park and the Carpathian Biosphere Reserve, the maximum permissible load on the 
Danube ecosystem should be the one when the WPI does not exceed the natural background value (long-term 
annual average standard), calculated in accordance with the dependences described above. These functional 
models for the monitoring assessment of the anthropogenic load on conservation areas should be used. 


To process the monitoring base of the second large Dniester hydro-ecosystem, the Composite Potential 
Quality Indicator (2) was used, and the close regression functional dependence between the hydro- 
ecosystem quality standard (long-term annual average) of the whole Dniester basin in the Carpathian region 
and the terrain altitude was obtained (Figure 6). 
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Figure 6: Functional dependence between the water quality standard of the Dniester ecosystem and 
terrain altitude 


The calculated CPQI is a relative value that depends on the level of anthropogenic load. Even though the 
Dniester ecosystem is constantly polluted at most observation points (even with positive CPQI values), the 


close dependence between the Composite Potential Quality Indicator and the altitude is evident and is 
described by the curvilinear regression equation: 


CPQI=1.095+8.18H*° (11) 
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In this case, Fstat=253.4 (calculated)>Fstat=3.94 (table value). Thus, the correlation between the 
properties is not accidental, it is significant. 


The functional dependence obtained for the Dniester ecosystem in the Ukrainian Carpathians by 
interpolation, 1.e., by means of the mathematical justification for the unknown values of the dynamic series 
of phenomena based on the established relationship between the CPQI standard and the altitude, was 
approximated by a continuous linear function (Figure 7): 


CPQI =1.002+2.705H’” (12) 


y=atbx"(2_5) 
r2=0.70220002 DF Adj r2=0.69641749 FitStdEn=0.4109969 Fstat=245.22769 
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Figure 7: Approximated functional dependence between the quality standard value of the water 
bodies of the Dniester ecosystem in the Carpathian region and the terrain altitude 


Thus, based on the acquired dependence, it is proposed to define the Composite Potential Quality Indicator 
standard value as a property of natural and technogenic safety at any point of the Dniester river ecosystem 
in the Carpathian region. 


Undoubtedly, the CPQI standard values at absolute altitudes below 300 m were influenced by the existing 
technogenic load on the ecosystem. The surface water quality monitoring point on the Tysmenytsia river is 
located in the town of Drohobych and 1s aimed to control the influence of sewage water from the towns of 
Truskavets and Boryslav. In Rozvadiv village, water quality of the Dniester ecosystem is monitored due to 
the impact of industrial wastewater from the town of Drohobych, and in Zhuravno village it was due to the 
impact of the sewage waters from the town of Khodoriv (Malovanyy et al., 2019). 


To determine the standard of the qualitative component of ecological safety of hydro-ecosystems of the 
Ukrainian Carpathians, it is proposed to use the Table 2. In the table, the norm of CPQI (for the natural state 
of hydro-ecosystems) could be estimated using the value of the height sample location. Up to an altitude of 
1,000 m above sea level, water quality can be determined at every 50 m distance, and at an altitude above 
1000 m it should be measured at every 100 m distance. 
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Table 2: Standard values of the Composite Potential Quality Indicator for different absolute heights of the 
hydro-ecosystems of the Dniester, Prut, Siret and Tysa rivers 
Hydro-ecosystem 


Absolute height, m 
Dnister 

O00 
150 
300 
450 
FO 2 
700 969 SB 


959 


Discussion 


Water quality assessment is a time-consuming process because it is based on the comparison of the average 
concentrations observed at the water quality monitoring points keeping in view the standards for each 
ingredient (Stasko and Buczynski, 2018). 


Following the principle of natural water unity, formulated by (Eigen and Schuster, 2012), the current system 
of water quality (drinking, waste, surface and underground waters) should be based on the classification 
built upon the indicators and the composition and properties (physical, chemical, biological) of water. It can 
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be collectively used to solve a wide range of problems associated with different types of management, 
utilization and protection of water resources. The system should meet simultaneously the ecological, 
hygienic and technological requirements (Zubaidah, Karnaningroem and Slamet, 2019). At present, there 
exists no such system. Hence, it was a complex multidisciplinary task (Water Framework Directive 
2000/60/EC’). 


Difficulties arise when it becomes inevitable to analyze the trend of water quality over several years, in 
different parts of a water body, or to compare the water quality of different water bodies contaminated by 
variety of pollutants, or to identify the trend of water quality over time (Strzelczyk and Steinhoff- 
Wrzesniewska, 2020). Thus, a need to develop a methodology for integrated water quality assessment arises 
(European Union, 2017). 


This study focusing two major hydro-ecosystems of the Carpathian region, the Dniester and Danube rivers, 
used a methodology of integrating natural water quality assessment with the functional dependencies 
between quality indicators and altitudes. It is established that the relation is logical: with the increasing 
altitude, the level of anthropogenic load decreases, and the quality of natural waters improves. It has been 
observed that the Ukrainian Carpathians are low mountains, quite populated. In small settlements, there 1s 
no centralized water supply and sewage treatment. As a tourist region, the Ivano-Frankivsk oblast has the 
population of 1.4 million and is visited by more than 2 million tourists annually. A big ski resort, “Bukovel’, 
is situated in the area (Kinash et al., 2019). Both the high-density population and the influx of tourists add 
to the pollution load of watercourses in the region. 


Based on the above findings, it can be inferred that the self-purification of natural waters eliminates diffused 
pollution in the waters flowing through mountainous part of the Ukrainian Carpathians. The comparison of 
method (2) with the method (1) leads to the following conclusions. The methodology (2) is more flexible, 
which evaluates the buffering capacity and allows the use of a simple prioritization method; that is, certain 
regions or areas of ecosystems that meet certain environmental quality standards can be considered as 
reference having no further anthropogenic pressure, while other areas of hydro-ecosystems can be ranked 
and evaluated using the mathematical symbols and values of CPQI. 


Finally, this method is proved to be effective in assessing the level of anthropogenic load on water bodies. 
The obtained dependencies can be used as territorial background standards for the ecological assessment of 
the basin ecosystems. This study proves to be important to justify and improve the Ukrainian water 
monitoring methodologies that are required to be adapted to EU laws. 


Conclusions 


For the first time, the patterns of changes in the ecosystems’ quality parameters established by means of 
developing functional dependences among the ecological standard values of the composite quality indicators 
of natural waters and the river length and terrain altitude. It is performed by using the statistical processing 
of the data of quality monitoring of the upper reaches of Dniester and Danube ecosystems in the Ukrainian 
Carpathians for the period from 2001 to 2019. The obtained dependencies can be used in the form of the 
territorial background standard values of ecological indicators of basin ecosystems. It will help substantiate 
scientifically the permissible levels of the anthropogenic pressure e.g., pollution load. 
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